Transport Properties near the z=2 Insulator-Superconductor Transition 
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We consider here the fluctuation conductivity near the point of the insulator-superconductor tran- 
sition in a system of Josephson junction arrays in the presence of particle-hole asymmetry and Ohmic 
dissipation. The transition is characterized by the dynamical critical exponent 2 = 2, opening the 
possibility of the perturbative renormalization-group (RG) treatment. The coupling to the Ohmic heat 
bath, giving the finite quasiparticle life-time, leads to the non-monotonic behavior of the dc conductivity 
as a function of temperature in the leading logarithmic approximation. 



A. Introduction 

In granular superconducting thin films and fabricated 
2D Josephson junction arrays (JJA), phase coherence is 
destroyed at zero temperature, whenever the intergrain 
Josephson tunneling of Cooper pairs becomes smaller 
than the Coulomb interaction between grains. [Q-^. 
The occurrence of this insulator-superconductor quan- 
tum phase transition (1ST) can be best seen from con- 
ductivity measurements. At the transition point, the re- 
sistivity is observed to be temperature-independent and 
close in value to the quantum of resistance for charge 
2e particles Rq — ctq^ = /i/4e^. There are many theo- 
retical works devoted to the explanation of this remark- 
able result |p[-pl|. Cha et.al., for instance, calculated the 
zero-temperature collisionless conductivity at the transi- 
tion point within the effective 3D Ginzburg-Landau (GL) 
action for an Af-component bosonic field that can 
be derived from the commensurate 2D Bose-Hubbard 
model. The influence of the quartic term in this ap- 
proach, considered up to the first order in 1/M, was 
shown to reduce the mean-field value of the conductivity 
(7r/8)crQ = ne^ /2h by 36%. Fazio and Zappala calculated 
this conductivity applying the dimensional regularization 
in e = 3 — c? p|. However, as was shown later by Damle 
and Sachdev |[11| , the collisionless zero-temperature con- 
ductivity is not the conductivity that is measured exper- 
imentally at a finite temperature. This conclusion stems 
from the non-commutativity of the limits w — > 0, T = 
and = 0,T — !■ in the general expression for a{Lu,T) 
at the transition point. The former limit describes colli- 
sionless transport, while it is the latter limit that must 
be taken to assess the experimentally-relevant collision- 
dominated conductivity. With collisions neglected, the 
finite temperature dc conductivity is singular on the in- 
sulating side, so the inclusion of quasiparticle damping 
is necessary for its regularization. Collisions are prop- 
erly accounted for utilizing the quantum kinetic equation 
which represents, in general, a complicated non-linear in- 
tegral equation for the distribution function of thermally 
excited quasiparticles. The kernel of this equation is de- 



termined by the quartic non-linearity in the GL action 
responsible for collisions between quasiparticles. The ap- 
proximate methods, such as 1/AI or e = 3 — d expansion, 
should be used to solve the kinetic equation close to the 
transition point Another way of regularization 

of the singular Drude conductivity is the inclusion of the 
phenomenological Ohmic dissipation. Experimentally, 
the latter plays a central role, because in combination 
with disorder it can lead to the temperature-independent 
resistivity as T ^ iQ. From the theoretical point of 
view Ohmic dissipation is always a relevant perturbation, 
affecting thus the functional form of universal quantities 
near the 1ST point PJT^]. The mentioned sources of dis- 
sipation lead to different temperature dependences of the 
conductivity, and, generally speaking, should be treated 
simultaneously p^ ]. 

Thus far, most of the approaches have used as a start- 
ing point the particle-hole symmetric GL action that is 
space-time isotropic at T = and belongs to the uni- 
versality class with the mean-field dynamical exponent 
z = 1. However, the behavior in the critical region can be 
considerably affected by a term that breaks the symmetry 
between particles and holes. The physical picture lead- 
ing to the modified action is realized in a system of JJA 
in the presence of uniformly-distributed frustrating offset 
charges Qx ^ 2e that cannot be eliminated by Cooper pair 
tunneling |^,|l6| . The microscopic quantum phase model 
describing this situation is also known to be equivalent 
to the clean incommensurate Bose-Hubbard model p7| . 
The 1ST in this case is characterized by the z = 2 dy- 
namical critical exponent. Thus, the quartic term in the 
GL action is marginally irrelevant, and there exists the 
range of parameters T and 8 (the deviation from the zero- 
temperature critical point) that permits all interesting 
physical quantities to be calculated in the perturbative 
renormalization-group (RG) treatment. p"^-|20[|. In this 
region the static quantities can be regarded in the lead- 
ing approximation as some universal functions of 8 and T, 
independent of the quartic interaction strength, u. This 
region persists as long as 8 and T are logarithmically or 
even double logarithmically small. 
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In this paper we use this approach to calculate the dc 
conductivity close to the point of the 1ST, governed by 
z — 2^ within leading logarithmic accuracy. It is shown, 
using the Kubo formula, that in the absence of any dis- 
sipative mechanism, the real part of finite temperature 
conductivity is singular and non-universal on the insulat- 
ing side. Its regularization by inclusion of the coupling 
to the Ohmic heat bath leads to the logarithmic temper- 
ature dependence in the quantum critical (QC) regime. 
The conductivity increases in the leading approximation 
as lnln(l/r) with T ^ right above the transition, and 
decreases monotonically in the quantum disordered (QD) 
regime. Consequently, we claim then that crossing over 
from the QD to the QC regime with decreasing temper- 
ature is accompanied by a somewhat re-entrant behavior 
of the dc conductivity. The regularization that results 
from the finite lifetime of quasiparticles, is possible only 
once the Umklapp scattering is included. This is a con- 
sequence of the Galilean invariance of the action, which 
explicitly describes the system of interacting charges of 
one sign. These processes ensure that the inverse life- 
time of quasiparticles, 1/tu ~ e~^^'^ , leads, in the ab- 
sence of other sources of dissipation, to an exponentially 
large fluctuation conductivity on the insulating side. We 
discuss the applicability of this approach to the descrip- 
tion of the finite-T charge transport near the 2D metal- 
superconductor transitions. 



B. Renormalization Group Analysis 

The general form of the GL functional that models the 
behavior near the 1ST point in the presence of uniformly- 
distributed offset charges is derived to be S 



d^r / dr 



V-H— A(r,r) ) V*(r,r) 



V - ^A(r, r) ) V(r, r) 



W |9.^(r,T)|V<5 l^(r,T)|' + ^ |V(r,T)|*} 



that in 2D at finite temperature, the conductivity obeys 
the scaling [p2|j23| 



+Fdi> 



(1) 



where A(r,T) is the vector potential, due to applied 
electric field, e* — 2e, and 5 is proportional to the in- 
verse correlation length. In Fourier space, the dissipation 
term, Fdis = '?Z]k,w„ kn||V'(k,t^n)P corresponds to the 
Ohmic model of Caldeira and Leggett |2^. Parameters 
AC and A (linear in the value of offset charges 1^) measure 
the strength of quantum fluctuations. We will regard 
k/A — 0(1). The term with k is clearly irrelevant in the 
region of interest, so we can neglect it and measure all 
quantities, having the dimensionality of energy, in units 
of A. The key idea behind the calculation of the general, 
frequency-dependent conductivity in the critical region is 



(2) 



In the momentum-shell RG, I* denotes the scale at which 
the effective size of the Kadanoff cell is on the order of 
the correlation length. We assume without loss of gen- 
erality that the scaling stops when 5{l*) — 1. The fre- 
quency in the above equation scales trivially, lo{1) = oje^\ 
although we will be interested only in the dc conductiv- 
ity, most commonly measured in experiments. The dis- 
sipation term with r], preserving the z = 2 universality 
class, should be also incorporated in the RG equations. 
However, because any non-zero rj profoundly changes the 
analytic structure of a quasiparticle propagator, it is con- 
venient to consider the cases of zero and non-zero rj sep- 
arately. 

77 = 0: The one-loop finite-T RG equations for this 
problem in 2D are of the form |T8| 



dTjl) 
dl 



dSjl) 
dl 



du{l) 
dl 



25{l) 



2K2u{l) 



exp[(l + 5{1))IT{1)] 



{2 - z)u{V) - C[5{l),T{l)]u{lf 



(3) 



(4) 



(5) 



where the coefficient 



cm).T{i)] 



2{l + 5{l)) 
5K2T{1) 
il + 5{l)) 



T{1) < 1 

2 m » 1 



(6) 



Here K2 — (l/27r), and in the one-loop approximation 
we can assume z = 2. The solution of the above system 
is qualitatively different in the quantum disordered and 
quantum critical regimes. 

Quantum disordered regime: We are in the QD regime, 
if the inequality 6 ^ T between the bare parameters 
holds. With sufficient accuracy the term with u{l) in the 
righthand side of Eq. ^ can be neglected, and we obtain 
in this regime, that /* = | In (1/(5), giving 



T 



Stt 



(7) 



The asterisk denotes, that the parameters under scaling 
must be taken aX I ~ I* . The result for u* is obtained 
with logarithmic accuracy, i.e we consider I* to be large. 
This condition ensures the validity of the perturbation 
theory. 

Quantum critical regime: This region is characterized 
by the opposite condition, 5 <^ T. The integration of 
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the RG equations in this case consists of two steps. In 
the first step, we integrate over I from to such that 
T{1) = 1. We have then 5(0 = (5/T + 0(l/ InT-i)) < 1. 
So, the scahng must continue into the second step in 
which T{1) 3> 1. In this, classical region, the strength 
of interactions is measured not by u{l), but by v{l) ~ 
u{l)T{l) U^, and for Z > / we must consider the equations 



dvjl) 
dl 



2v{l) 



{l + 5{l)y 



dl 



l + 5{l) 



(8) 



(9) 



Assuming, that throughout the scaling v{l) <C 1, that is 
we are working in the weakly classical region, we obtain, 
using v{l) = 47r/f, that v{l) « (47r/f)e2('"'l Substitution 
of this result into Eq. (|9|) and subsequent integration 
gives us with the logarithmic accuracy the equation for 
I* in the QC regime: 



5{r 
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l)e 



2(1* -I) 



1. 



(10) 



Solving this within double logarithmic accuracy, we find 
that 



I* 



1 



■In 



^lnln(i; 



T* 
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(12) 



The validity of the above solution requires not only the 
condition ln(l/T) > 1, but also lnln(l/T) > 1. 

ry ^ 0: The presence of the term ri\ujn\ in the RG 
equation for 5 leads to the divergency in the sum over 
frequencies. So the momentum-shell RG with the upper 
frequency cutoff F should be used , and the obtained 
one-loop equations are of the form 



ddjl) 
dl 



^26il) + f<^'\d{l),Til))u{l) 



du(l) 
dl 



~f^^\S{l),T{l))u{l)' 



(13) 



(14) 



where /(^' and /("^^ are some rather robust functions of 
S{1), T{1) and 77, that we don't write out here. It is es- 
sential, that the results for T* and /* remain the same in 
both the QC and QD regimes, if we are interested only 
in the leading logarithmic approximation. However, the 
subdominant temperature dependent corrections appear 
to be different, u* and v* are also logarithmically small, 
though have the cutoff dependent numerical prefactors. 



C. Transport Properties 

We see, that upon scaling, for small enough T and (5, 
the parameters u* and v* remain small, and in the first 
approximation one can calculate the conductivity with 
the help of the Kubo formula. 



a'Q/3(jWn) = 



dT- 



52 \nZ 



applied to the Gaussian part of Eq.(^. The standard 
calculations for the longitudinal conductivity, which we 
denote simply as ct, lead to the result 



2fc2G(k,0 



(G(k, uj*^) - G(k, uj^^ + uj*)) , 



(15) 



where G(k, lj*) = (iw* + el)^^ is the usual Matsubara 
Green function. The rescaled temperature T* and the 
energy of quasiparticles e^. = 1 + are employed in 
the righthand side in accordance with Eq. (||) {uj^ — 
2TnnT*). We must perform then the analytical continu- 
ation to real frequencies after doing the summation over 
frequencies to^. 

If the absence of any dissipative mechanism, we obtain: 



(11) a{io*) = 



— +7rS{Lu*) 
uo* 



2hT* 



k^dk 



sinh 



l + fc^ 
2T* 



(16) 



Though the integral over k is calculable exactly, we state 
here only the results in the limiting cases of the QD and 
QC regimes: 



(7(lu) = h tt5(lu) 



T 



p-S/T 

In In In ^ 



QD 
QC 



(17) 



The expression in the QC regime is written with triple 
logarithmic accuracy. We see that on the insulating side, 
the real part of the coUisionless conductivity does not 
have any regular contribution, as — > 0. The singu- 
lar part disappears as well at T = 0. Moreover, the 
form of the frequency dependence in Eq. (|lj) is charac- 
teristic for the conductivity of a superconducting phase. 
This unphysical result suggests that the finite life-time 
of quasiparticles should be accounted for in the correct 
determination of transport properties. 

The inclusion of the non-zero Ohmic dissipation pro- 
vided by the last term in Eq. (|l|), changes the analytical 
properties of Matsubara sums, and the analytical con- 
tinuation in Eq. ( p^ ) should be performed by contour 
integration 



huj* 



POO POO 

/ k'dk / coth —dz [(G«(z) - G^(z)) 

>/ tJ — 00 



[G«(z) + G^{z) - G^(z + UJ*) - G^{z - LU*)]] , (18) 
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where the retarded and advanced Green functions 
G^^"^{z) — (z + =F ir]z)~^ have been introduced. In 
the hmit w — > 0, expansion of the corresponding Green 
functions and subsequent integration by parts yields the 
dc conductivity: 



k^dk 



dx 



"^T^h Jq J-oo sinh^ x 

8iTjT*x)^ 

[(e* + 2xT*)^ + {2rjT*xf 



(19) 



In the QD regime T* — T / 6 <^ 1, and, as long as 
77 <C the integral over x is determined by the two 
competing contributions: from the region near the min- 
imum of denominator at xq ~ — e^/2T* and from the 
vicinity of a; = 0. Evaluating those contributions and in- 
tegrating subsequently over momentum fc, we receive the 
total conductivity: 



2e2 
h 



T 
Sr] 



-S/T 



27r f-qT 



(20) 



As can be seen, the first term is dominant for very small 
values of dissipation, 77 < {S/Ty^^e~^^^'^ , while for larger 
rj only the second contribution should be retained. In any 
of that cases the conductivity monotonically decreases as 
T ^ 0. 

In the QC regime the double logarithmic accuracy im- 
plies that T* = lnln(l/T) 3> 1. With the same accuracy, 
the main contribution to the integral over x comes from 
X <^ 1, allowing us to set x^/sinh^x w 1. We obtain 
then after simple integrations, assuming rj < I, that 



— ^ In In — . 

; T 



(21) 



This suggests, that the conductivity is not universal in 
the QC regime and increases upon lowering the temper- 
ature, albeit very slowly, as a double logarithm of 1/T. 
One can conclude that, because of the crossover from the 
QC to the QD region upon lowering temperature for any 
non-zero 6, it may be possible to observe a non-monotonic 
behavior of the resistivity as a function of temperature 
with a dip when T ~ 5. Though the derived results are 
strictly applicable in rather tiny regimes, where S and T 
are logarithmically small, the results found here may be 
observable beyond those boundaries. 

Now consider the role of mutual scattering of quasipar- 
ticles as the other possible source of dissipation. Eq.(|^) 
suggests that our action is Galilean invariant, describing 
the system of charge carries of one definite sign. In the 
Hamiltonian formalism, our system is identical to that 
of weakly-interacting Bose particles in 2D j24j . The lin- 
earization of the corresponding collision integral in the 
kinetic equation in a small correction to the equilibrium 
distribution, proportional to (kE), 



m-^)fiKt)^i[fik,t)], (22) 



d_ 

dt 

reveals, that /[/(k, t)] turns to zero as a consequence of 
momentum conservation. This situation is different from 
the z = 1 case, in which the system consists of colliding 
particles and holes having the same energy |ll[. This 
does not mean, however, that the mutual collisions of 
quasiparticles in the z = 2 case can not lead to the finite 
lifetime. The appropriate way to see this is to go beyond 
the continuum limit and recall that the real experimental 
systems, such as JJA, have a lattice periodicity. This 
suggests that in case of a square lattice, the quasiparticle 
energies Ck — 4 — 2 (cos 

~1" cos ky^ are periodic functions 
of the quasimomentum k, and the Umklapp processes 
need to be taken into account. The corresponding inverse 
scattering time 1/tu can be estimated from 

— ocu"^ n(eki)(l + "-(ek2))'5(ek + Cki - Cka - ekg) 

TU J 

d^ki d'^k2 d'^ks 



•(27r)^(5(k + ki-k2-k3-b) 



(27r)2 (27r)2 (27r)2 



(23) 



where the integration runs over the first Brillouin zone, 
and b is the reciprocal lattice vector. The non-zero b im- 
plies that the conditions imposed by 5-functions on quasi- 
momenta and energies of colliding quasiparticles can both 
be satisfied only if one of the incoming and one of the out- 
going quasimomenta come from the interior of the first 
Brillouin zone. Because the main contribution, as can 
be seen from the Kubo formula, comes from small fc, ^2 
must be large. Physically, this means that the damping 
of the critical fluctuations of the order parameter arises 
from the collisions with the non-critical modes. One can 
easily estimate that I/tu oc e'^^"^ where the prefactor 
varies slowly S and T, and A is some numerical constant 
determined by the concrete form of the first Brillouin 
zone and the quasiparticle dispersion law. This inverse 
scattering rate, valid in both QC and QD regimes, is very 
small near the quantum phase transition point. This im- 
plies that in the absence of any Ohmic dissipation rj the 
conductivity a oc T(j oc e^^"^ , and hence exponentially in- 
creases with temperature. The latter conclusion can be 
verified by making the substitution rjxT {\/2tu) in 
Eq. ( p^ ) . Such behavior suggests that Ohmic dissipation 
plays the dominant role near the 2D 1ST point in the 
presence of charge frustration and reinstates the insulat- 
ing behavior. For JJA, fabricated from s-wave supercon- 
ductors, weak non-magnetic disorder is non-pair break- 
ing. Hence, as long as the pairs remain intact, it can not 
lead to dissipative terms in our action |^^. Ohmic dissi- 
pation arises if the superconducting grains are embedded 
in a conducting environment ||2^ in which case Eqs. ( |20| ) 
and (|2^) give the conductivity due to fluctuations of the 
superconducting order parameter. 

Conventional disorder leads, however, to dissipation 
in a d-wave superconductor because in this case it is a 
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pair-breaking perturbation |2^]. The finite temperature [27] I.F. Herbut, Phys. Rev. Lett. 85, 1532, (2000). 
fluctuation conductivity near this disorder-driven metal- 
superconductor quantum phase transition can be investi- 
gated using the same action (||) , provided we are not too 
close to the transition point where the disorder affects 
the critical properties. Depending on parameters of the 
microscopic Hamiltonian, rj can be both small and large, 
compared to unity |^^. If 77 < 1, the conductivity in 
the QC regime is given by Eq. and we should use 

the last term of Eq. ( [20| ) for the QD regime. If 77 ^ 1, 
the relevant conductivity was discussed using N — 00 
approach Q. 
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